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Abstract. The main point of this communication is that there is a small non-negligible 
amount of eddies-outliers/very strong events (comprising a significant subset of the 
tails of the PDF of velocity increments in the nominally-defined inertial range) for 
which viscosity /dissipation is of utmost importance at whatever high Reynolds number. 
These events contribute significantly to the values of higher-order structure functions 
and their anomalous scaling. Thus the anomalous scaling is not an attribute of the 
conventionally-defined inertial range, and the latter is not a well-defined concept. The 
claim above is supported by an analysis of high-Reynolds-number flows in which among 
other things it was possible to evaluate the instantaneous rate of energy dissipation. 

There are a variety of models of higher statistics that have meager or nonexistent de- 
ductive support from the NS equations hut can he made to give good fits to experimental 
measurement^. These include 'explanations' of what is called anomalous scaling ob- 
served experimentally for higher-order structure functions of velocity and temperature 
increments, such that their scaling exponents = p/3 — < p/3 are nonlinear concave 
functions of the order p. Starting with refined similarity hypotheses by Kolmogorov'^ and 
Oboukhov'^, numerous phenomenological models have been proposed to describe these de- 
viations considered as the major manifestation of intermittency in the inertial rangd^-Hl. 
The dominant of these models has been the multi-fractal formalism^, others claimed the 
Reynolds number dependence as responsibl^^H^^. The common in all these approaches is 
the basic, widely accepted premise that in the inertial range, the viscosity plays in prin- 
ciple no rol^^ so that nonlinear dependence of the algehraic scaling exponents (p on the 
moment order p is a manifestation of the inertial-range intermittenc^ with the inertial 
range defined as 77 <^ r <^ L (with rj being the Kolmogorov and L — some integral scale). 
Thus the issue is directly related to what is called inertial (sub)range and how inertial it 
is. 

The main point of this communication is that there is a small non-negligible amount 
of eddies-outliers/very strong events (comprising a significant subset of the tails of the 
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FIG. 1. (a) Scaling exponents of structure functions at Rex ~ 10^ for the longitudinal velocity 
component corresponding to the full data and the same data in which the strong dissipative 
events (when at least at one point x or x + r the instantaneous dissipation e > q{e)) with various 
thresholds q were removed, (b) Scaling exponents for the strong events themselves. 

PDF of Aui{r) in the nominally defined inertial range 77 <^ r <^ L) for which viscos- 
ity/dissipation is of utmost importance at whatever high Reynolds number. In other 
words, the inertial range is ill-defined in the sense that not all, but almost all statistics of 
Aui{r) is independent of viscosity. As long as one deals with low-order statistics of Aui{r) 
(as Kolmogorov did) this is of little (but not always negligible) importance. However, it 
appears that these events contribute significantly to the higher-order structure functions 
and thereby a non-negligible contribution to the higher-order structure functions is dom- 
inated by viscosity. In other words, the 'anomalous scaling' as exhibited by the behavior 
of higher-order structure functions is to a large extent due to significant contribution of 
viscosity/dissipation in the inertial range as commonly defined. The higher the order of 
the structure function, the stronger is the contribution due to viscosity (i.e. from the 
tails of the PDFs of Aui{r)) and the weaker is the 'inertial' contribution (i.e. from the 
core of those PDFs) to the structure function. Thus it seems problematic to speak about 
inertial-range behavior of higher-order structure functions. 

The support for the above view comes from a recent analysis of high- Reynolds-number 
data in field experiment The experimental facilities and related matters are de- 
scribed in these papers and references therein. We give here a very brief reminding on 
these. 

The measurement system, developed by the group of Prof. Tsinober, consists of the 
multi-hot-wire probe connected to the anemometer channels, signal normalization device 
(sample-and-hold modules and anti-aliasing filters), data acquisition and calibration unit. 
The probe is built of five similar arrays. Each calibrated array allows to obtain three 
velocity components "at a point". The differences between the properly chosen arrays 
give the tensor of the spacial velocity derivatives (without invoking of Taylor hypothesis), 
temporal derivatives can be obtained from the differences between the sequential samples. 



2 



10000 



5000 




100 



200 



300 



400 



500 



600 




200 300 400 

Sample number, from position 48000 



600 



FIG. 2. Example of simultaneous time series of the squared magnitude of the rate of strain 
tensor, s^, proportional to the dissipation e (top) and the velocity increments, Aui = ni(x + r) — 
ui{x) for r = 400r/ (bottom). The marked segments correspond to the strong events, selected 
with the value of the threshold q = 12 equivalent to the value of ~ 4, 000. It should be noted 
that the first two marked segments are considered strong events because the value of reaches 
the threshold at the point x + r. 




FIG. 3. PDFs of the increments of the longitudinal velocity component for the same data as 
in Fig. [T] (a) r/r/ = 40 corresponds to the lower edge of the inertial range, (b) r/r] = 400 is 
deep in the inertial range. 



The calibration unit produces a jet with variable velocity magnitude and variable 
angles around two orthogonal axes. The probe is located in the jet core, the values of 
velocity magnitude and angles are recorded together with the readings of the anemometer 
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channels and later approximated by polynomials, used for obtaining velocity components 
from the measured hot-wires data. 

A selection of results is shown in Fig. [T] — Fig. [6j 

Fig. [l (a) shows the scaling exponents of structure functions Sp{r) up to order 8 
corresponding to the full data and the same data in which the strong dissipative events 
with various thresholds were removed. By an event we mean here a velocity increment, 
Aui = Ui{x + r) — ui{x). It is qualified as a strong dissipative event if at least at one 
of its ends {x,x + r) the instantaneous dissipation e > g(e) for g > 1. We have chosen 
q = 3,6,12 and 20. This corresponds to the instantaneous Kolmogorov-like scales 0.76, 
0.64, 0.54 and 0.47 of the conventional Kolmogorov scale rj based on the mean dissipation 
(e). It is seen that the removal of the strong dissipative events results in an increase 
of the exponents (p. For example, with the removal of the dissipative events between 
the threshold 3(e) (0.767;) and 6(e) (0.647;) the dependence of (p on p becomes pretty 
close to the Kolmogorov p/3. The strong events/outliers themselves have different scaling 
properties (Fig. [l](b)). The time series in Fig. [2] illustrate the selection of strong events. 

The next example in Fig. [3] shows that indeed the removal of the strong dissipative 
events results in narrowing of the tails in the PDFs of Aui{r). 

As an additional illustration we show in Fig. |4] two examples of histograms of the in- 
crements of Ui for the whole field, with removed strong dissipative events for the threshold 
q = 3 and the dissipative events themselves for the same threshold, for the same data as 
in Fig. [TJ 

The effect of the removal of the strong dissipative events is obviously much stronger 
for higher-order structure functions. For example, there are only 5% of dissipative events 
(Fig. |5] (a)) for g = 6 sitting mostly at tails of the PDF of Aui{r) for r /rj = 400 (i.e. 
deep in the 'inertial' range), which contribute about 36% to the total dissipation (Fig. Is] 
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FIG. 5. Percent of the strong dissipative events as defined in the text (a) and their contribution 
to the total dissipation (b) as a function of the threshold q for various separations r. 
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FIG. 6. Contributions of the strong dissipative events, as defined in the text, to the eighth-order 
structure function (a), the second-order structure function (b) and the third-order structure 
function (c) as a function of the threshold q for various separations r. 



(b)). These events contribute nearly 60% to the value of 5'g(r) at Rex ~ 10"^ (Fig. |6](a)). 
These same events change the 5'2(r) by about 11.5% (Fig. [g] (b)), but contribute about 
9% to S'^{r) (see Fig. |6|(c)). 

It is noteworthy that the data used herd ^^ l ^'^ l was somewhat spatially underresolved, 
1 ^ St]. This means that the conclusions are to some extent qualitative. However, with 
properly resolved data the strong dissipative events, lost in the underresolved ones, would 
somewhat enhance the tendencies just described above. This is in agreement with the fact 
that essentially the same results are obtained using the same data smoothed over up to 
eight sequential samples. Additional support comes from referenc^, indicating that the 
underresolved data reproduce faithfully the flow at scales about two times smaller than 
those resolved (~ O.Gr]) at least as concerns the instantaneous dissipation rate. Finally, 
using enstrophy u"^ and/or the surrogate {dui/ dxiY as a criterion for the threshold instead 
of the true dissipation e gives the same qualitative (but not quantitative) results. 
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Summarizing, the main point is the distinction between roughly two kinds of events 
(in the nominal inertial range r] r <^ L), contributing to the value of the velocity 
increment. One is represented by the core of the PDF of Awj, and the other — by the 
outliers/extreme events (comprising a significant subset of the tails of the PDF of Awj). 
They have not only different statistical properties (such as scaling if such exists), but 
also are of different nature in the sense that the former exhibit 'inertial' behavior as 
reflected in the slopes of low-order structure functions, whereas the latter are dominated 
by viscous effects as seen in the slopes of higher-order structure functions. Removal 
of these highly-dissipative events brings the dependence of Cp on P pretty close to the 
Kolmogorov p/3. Thus the anomalous scaling is not the attribute of the inertial range. 
Our results leave little doubt that the strong dissipative events contribute significantly to 
the anomalous scaling of higher-order structure functions. However, there are other effects 
which are expected to contribute to 'anomalous scaling' such as a variety of nonlocal effects 
understood in a broad sense as direct and bidirectional coupling/interaction between large 
and small scale^lEH The quality of our data does not allow to address properly this and 
similar issues. This is a matter of far more precise and well-controlled experiments which 
among other things require information at high Reynolds numbers with sub- Kolmogorov 
resolution. 

Along with the fact that velocity increments (let alone structure functions and their 
scaling if such exists) are not the only objects of interest and do not constitute a repre- 
sentation basis for a flovJ^, they are not a good object to define a perfect inertial range. 
Such a definition seems to be not possible in principle due to a variety of nonlocal effects 
as mentioned abov^3iH 

A special remark is about the contribution of the dissipative events as defined/described 
above to the 4/5 law. These events do contribute to the 4/5 law, and removing them 
leads to an increase of the scaling exponent above unity, see Fig. |l|(a) and Fig. |6|(c). An 
important point here is that the neglected viscous term in the Karman-Howarth equation 
does not contain all the viscous contributions. Those which are present in the structure 
function 5*3 itself remain and keep the 4/5 law precise. It this sense this law is not a pure 
inertial law. In fact, the contribution of the strongly-dissipative events is non- negligible 
also in the core of the PDFs of Ami (but not dominating as in their tails) as can be seen 
from Fig. |4| 

Among the main challenges for future work is an experiment similar to that described 
in reference£3Ill but with sub-Kolmogorov resolution. This includes also the issue of 
passive scalar. So far, we have pretty crude qualitative results (due to poor resolution 
and quality of the data) concerning the passive scalai'^, which show the same trends as 
described above and which raise similar questions concerning the anomalous scaling of 
passive scalars. 
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